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by 

Joseph C. Cain and Shirley J. Cain 
Goddard Space Flight Center 


INTRODUCTION 

This paper summarizes some of the computations that were made at the International Association 
of Geomagnetism and Aeronomy (lAGA) Symposium in Washington, D. C., October 22-25, 1968, 
that led to the resolution by the Working Group on the Analysis of the Geomagnetic Field (Reporter, 
A. J. Zmuda) to propose an International Geomagnetic Reference Field (IGRF). The basic require- 
ments established by Dr. Zmuda following the discussion at previous meetings called for the IGRF to 
consist of no more than 80 spherical harmonic coefficients of internal origin, epoch 1965.0, each be- 
ing tabulated together with its first time derivative. These coefficients were to represent true spherical 
harmonics describing the field, not “quasi-spherical” harmonics resulting from derivations neglecting 
the oblateness of the earth. Further, only sets of coefficients submitted to the Working Group on or 
prior to March 15, 1968, were to be considered. 

These sets of spherical harmonic coefficients are given in Table 1 . They are each updated to 
1965.0 and are limited to an n* (maximum degree n and order m of the spherical harmonic coeffi- 
cients) of eight. Of the sets given, all except those in Tables 1 (g) and 1 (h) take into account the ob- 
lateness of the earth in their derivation. Most of the field descriptions appear in the World Mag- 
netic Survey (WMS) Volume (Zmuda, 1971). However, a few have also been published separately, 
as follows. 


Table 

Field Model 

Reference 

1(a) 

GSFC( 12/66) 

Cain et al., 1967 

Kg) 

USC&GS* 

Hurwitz et al., 1966 

1(h) 

RGO-1 (LME)* 

Leaton et al., 1965 


^United SUtes Coast A Geodetic Survey. 

^Royal Greenwich Obaervatocy (RGO) model 1, tMsed on Leaton, Malm, and Evans (1965). 
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TABLE 1(c) TABLE 1(f) 
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TABLE Kg) TABLE 1(h) 


TEST DATA 



Although no explicit formula was agreed upon prior to the meeting for the derivation of an 
IGRF, there was an understanding that the model had to correspond to the available survey data. 

Since the epoch of this IGRF was 1965, data were arbitrarily cut off at 1961, a year chosen so 
the results would not be too heavily weighted by observations prior to 1965. Testing was done on all > 

data available since that date. These were divided into the major categories below. 

(1 ) Observatory annual means of surface magnetic fields, 1961-1967. 

(2) Surface magnetic surveys. This category includes land surveys, repeat stations, shipboard 
and ship-towed observations. 

(3) Aeromagnetic survey of Japan, 1965 (Nagata, 1966). 

(4) Aeromagnetic survey of Canada, 1961-1963. 

(5) Aeromagnetic survey of Scandinavia, 1965 (Eleman et al., 1969). 

(6) Project MAGNET worldwide (principally oceanic) airborne survey, 1961-1966 (USNCX), ^ 

1965). I 

(7) OGO 2 data, as available during magnetically quiet intervals, October 1965 to September 

1967. 

(8) OGO 4 data during magnetically quiet intervals from July to December 1967. 

(9) 1964-83c observations, 1964-1965 (Zmuda et al., 1968). 

(10) Cosmos 49 observations, 1964.8 (Dolginov et al., 1967). 

(11) Other airborne (towed proton-magnetometer) data. 

All of the nonsatellite data were obtained from the file prepared by the Geomagnetic Division of i, 

the U.S. Coast and CJeodetic Survey (E. Fabiano and S. Cain, WMS Volume). This file con- 
tained the coulributions from many separate organizations and survey groups and is constantly up- 
dated as new observations are submitted. This file was edited by rejecting those observations deviating 
by more than lOODy from the GSFC( 12/66) model (using n* = 10). This procedure was used to elimi- 
nate the highly anomalous data beyond about five times the root-mean-square (rms) deviation. Since 
all models were truncated to n* = 8 for testing, no particular advantage was given to GSFC( 1 2/66). 

This model was used since it fitted the data set best; hence, it requires the least elimination of data. 

The amount rejected was small, as seen in Table 2. 


, The OGO 2 and OGO 4 data (sampled every 30 seconds or at a spacing of approximately 200 km) 

I were initially selected from periods of time for which - 0. They were then fit with a special model 

I listed in Table 3(a) [POGO( 10/68)] employing 143 internal coefficients and their first time derivatives. 

I The distribution of deviations of the data from this fit was as follows. 


6 


|AF|7 0 10 20 30 40 50 60 70 100 200 600 Total 

Obs. 27,646 4218 589 141 23 26 6 2 9 4 32.664 

Since the distribution indicated that they were probably anomalous, the 15 observations over 7O7 
were rejected and the resulting rms deviation computed to be ly. The remaining 32,649 observations 
were included in the testing. 

The Cosmos 49 data were similarly treated by fitting with a special function and eliminating those 
data that deviated significantly from the rest. The data were prepared by the U. S. Coast and Geodetic 
Survey from the catalog (Dolginov et al., 1967) published by the Institute of Terrestrial Magnetism 
and Radiowave Propagation (IZMIRAN). These were sorted into time order and each fourth observa- 
tion fit with a series of 99 spherical harmonic coefficients by a model labelled COSMOS(9/68), listed 
in Table 3(b). Data deviating more than IOO7 from the fitting surface were rejected in the coefficient 
determination. The distribution of residuals from this model, COSMOS(9/68), is as follows. 

IAFI7 0 10 20 30 40 50 60 70 80 90 100 Total 

Obs. 1853 1243 648 271 93 41 23 18 19 15 138 4362 

The use of every fourth observation in the fit is adequate since each orbit then contains about 10 
observations for the shortest wavelength of the fitting function used {n* = 9 corresponds to 360*/9 = 
40®). Since the rms deviation of these data from the COSMOS(9/68) field was 2I7, the selection used 
for model testing consisted of those deviating by less than 6O7, a total of 16,554 observations from 
the approximately 1 8,000 originally available. 

The 1964-83c observations entered the testing unedited except for the rejection of one spurious 
point that gave a |AF| > IOOO7. 


Table 2-Nonsatellite data eliminated for AC> IOOO7. 


Data Type 

Component 

Observations* 

Data Rejected 

Number 

Percent 

(1) Observatory 

1984 

34 

1.7 

(2) Surface 

22,425 

204 

.9 

(3) Japanese Air 

1461 

6 

.4 

(4) Canadian Air 

9470 

27 

.3 

(5) Scandinavian Air 

6973 

1 


(6) Project MAGNET 

104,228 

401 

.4 

(7) Other Air 

1763 

9 

.5 


*ln tlm and enniing discussion a value of D, I, H, Z. or F is counted as one observation 
even though other values nuy have been measured at the same time and location. 
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Table 3(a)-POGO( 10/68) spherical harmonic coefficients and their time derivatives (Epoch 1960.0). 


n 

m 

X 

h 

X 

h 

n 

m 

X 

h 

X 

h 

1 

0 

- 304 ft 5.0 


25.42 


9 

0 

11.0 


- 0.24 


1 

1 

- 2163.3 

5791.0 

9.88 

- 4.66 

9 

1 

6.6 

- 20.4 

0.30 

- 0.38 

•> 

0 

- 1541.4 


- 23.90 


9 

* 

1.8 

14.4 

0.04 

0.16 

2 

1 

2976.3 

- 1977.2 

3.50 

- 7.07 

9 

3 

- 12.5 

0.6 

0.04 

0.67 

■y 


1607.5 

156.6 

- 2.14 

- 10.70 

9 

4 

15.8 

- 1.5 

- 0.40 

- 0.14 

3 

0 

1325.8 


- 5.59 


9 

5 

1.7 

1.4 

- 0.28 

- 0.83 

3 

1 

- 1983.7 

- 445.3 

- 11.52 

8.48 

9 

6 

2.6 

3.4 

- 0.57 

1.30 

3 


1316.9 

233.4 

- 4.41 

0 . ft 8 

9 

7 

8.7 

14.8 

- 1.32 

- 0.33 

3 

3 

842.0 

- 94.9 

2.87 

- 14.89 

9 

8 

5.1 

2.4 

- 0.14 

- 0.38 

4 

0 

959.1 


- 0.62 


9 

9 

- 2.4 

- 0.9 

0.50 

0.99 

4 

I 

819.6 

135.4 

- 2.51 

3.45 

10 

0 

- 2.6 


- 0.01 


1 

2 

486.4 

- 266.7 

- 1.22 

- 0.39 

10 

1 

- 2.0 

1.1 

- 0.09 

0.21 

4 

3 

- 372.4 

20.7 

- 2.95 

- 0.87 

10 

2 

1.0 

0.9 

0.12 

0.05 

4 

4 

256.2 

- 241.5 

0.86 

- 6.52 

10 

3 

- 5.5 

- 0.3 

0.19 

0.54 

5 

0 

- 234.3 


2.72 


10 

4 

- 0.7 

7.5 

- 0.20 

- 0 . 2 ft 

5 

1 

357.7 

16.9 

0.48 

0.05 

10 

5 

7.5 

- 2.3 

- 0.02 

- 0.30 

5 

■> 

233.9 

113.3 

3.17 

3.00 

10 

6 

7.8 

1.4 

- 0.43 

- 0.03 

S 

3 

- 21.0 

- 128.7 

- 2.46 

0.32 

10 

7 

1.7 

- 0.5 

- 0.33 

- 0.39 

5 

4 

- 147.1 

- 115.1 

- 0.89 

3.11 

10 

8 

- 5.3 

4.3 

1.03 

- 0.02 

5 

5 

- 45.2 

130.3 

- 3.15 

- 6.35 

10 

9 

1.3 

8.0 

0.22 

- 1.04 

6 

0 

49.1 


- 0.61 


10 

10 

- 2.7 

- 13.7 

0.46 

0.79 

ft 

1 

54.5 

- 9.6 

1.06 

- 0.26 

11 

0 

2.3 


0.03 


ft 

2 

4.8 

106.4 

0.62 

- 0.48 

11 

1 

- 1.8 

- 0.8 

0.11 

0.35 

6 

3 

- 249.1 

56.8 

3.95 

2.58 

11 

2 

- 2.1 

4.4 

0.05 

- 0.26 

ft 

4 

1.7 

- 27.2 

- 0.94 

- 0.80 

11 

3 

5.5 

- 0.1 

- 0.30 

- 0.17 

ft 

5 

- 3.7 

- 14.9 

1.49 

0.50 

11 

4 

- 1.5 

- 3.9 

0.01 

0.17 

ft 

6 

- 91.6 

- 4.3 

- 1.67 

0.82 

11 

5 

2.4 

- 0.6 

- 0.34 

0.18 

7 

0 

75.9 


- 0.89 


11 

6 

- 3.5 

1.8 

0.48 

- 0.50 

7 

1 

- 52.4 

- 57.9 

- 0.21 

- 0.87 

11 

7 

- 1.3 

- 3.2 

0.52 

0.23 

7 

2 

8.0 

- 25.0 

- 1.08 

- 0.46 

11 

8 

2.5 

0.8 

- 0.19 

- 0.34 

7 

3 

10.0 

- 0.8 

0.70 

- 1.02 

11 

9 

- 1.2 

- 5.9 

- 0.03 

0.37 

7 

4 

- 36.7 

6.3 

1.07 

0.25 

11 

10 

12.7 

- 1.7 

- 1.65 

0.22 

7 

5 

- 8.3 

9.5 

0.96 

1.88 

11 

11 

5.0 

10.5 

- 0.40 

- 1.55 

7 

6 

6.6 

- 11.7 

1.01 

- 2.43 







7 

7 

- 22.7 

- 37.6 

5.23 

2.32 







8 

0 

7.4 


0.61 








8 

1 

6.0 

10.1 

- 0.12 

- 0.15 







8 


- 8.1 

- 13.0 

0.87 

- 0.10 



m 


1960 ) 


8 

3 

- 9.2 

11.5 

- 0.33 

- 1.22 







8 

4 

- 0.8 

- 16.4 

- 0.14 

- 0.26 



h ( t ) 

• h * ii(t - 

1960 ) 


8 

5 

9.1 

5.5 

- 0.85 

0.15 







8 

ft 

- 11.4 

22.3 

0.99 

- 0.37 







8 

7 

7.9 

- 4.9 

0.81 

0.29 







8 

8 

35.1 

- 26.2 

- 4.98 

0.91 








TEST RESULTS 


Ihe various models were tested against the data sets both with the limitation of 80 coefficients 
and also using all coefficients if more were available. Table 4(a) illustrates for the GSFC(1 2/66) model 
the distribution of residuals using the. first 80 coefficients as well as the full number. Since the surface 
data were edited with this model using a ICOOy criterion* there can be no residuals above this figure 
with 1 20 coefficients. The effect of the truncation is to increase the rms residuals by 107*207 inde* 
pendent of their magnitude. Using 80 terms has only a small percentage effect on the surface data 
since magnetic anomalies account for a great deal of the scatter. The consequences for the satellite 
data arc more obvious as seen in the (XX) 2 results. Here the effect is to increase the number of resid* 
uals in the S 07-1007 range from S to 10 percent of the total data* and to push the number over lOOy 
from 1 to 3 percent. 
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Table 3(b)-COSMOS(9/68) spherical harmonic coefficients. 


n 

m 

g 

h 

n 

m 

g 

h 

1 

0 

-30415.2 


7 

0 

30.3 


1 

I 

-2143.1 

5721.6 

7 

1 

-52.4 

-70.2 

2 

0 

-1640.5 


7 

*» 

4.5 

-28.2 


I 

3001.9 

-2014.4 

7 

3 

5.4 

-14.2 


”» 

1556.8 

189.2 

7 

4 

-20.3 

13.3 

3 

0 

1211.1 


7 

5 

-9.1 

29.6 

3 

1 

-2033.8 

-388.6 

7 

6 

11.8 

-15.6 

3 

2 

1286.4 

258.0 

7 

7 

-12.1 

-15.4 

3 

3 

780.0 

-233.8 

8 

0 

15.5 


4 

0 

969.6 


8 

1 

11.1 

3.8 

4 

1 

816.3 

137.6 

8 

2 

-8.5 

-21.1 

4 

2 

487.2 

-301.9 

8 

3 

-9.3 

2.2 

4 

3 

-3867 

0.5 

8 

4 

-6.8 

-12.2 

4 

4 

253.6 

-186.4 

8 

5 

17.0 

-1.6 

5 

0 

-299.4 


8 

6 

7.7 

27.1 

5 

1 

348.5 

-0.5 

8 

7 

16.6 

-9.3 

5 


264.3 

106.7 

8 

8 

7.9 

-6.5 

5 

3 

-12.4 

-98.3 

9 

0 

-9.1 


5 

4 

-172.3 

-108.3 

9 

1 

7.2 

-29.3 

5 

5 

-35.5 

57.0 

9 

2 

II.O 

5.6 

6 

0 

57.6 


9 

3 

-14.6 

13.5 

6 

1 

69.2 

-19.1 

9 

4 

9.5 

-2.3 

6 

2 

5.9 

110.4 

9 

5 

2.5 

-5.1 

6 

3 

-228.1 

74.9 

9 

6 

0.7 

6.4 

6 

4 

6.4 

-56.0 

9 

7 

3.9 

9.6 

6 

5 

-19.4 

3.3 

9 

8 

3.3 

-1.7 

6 

6 

-158.9 

-33.8 

9 

9 

-2.3 

0.5 


These distributions were also calculated for each of the other test models, and the rms values 
compiled in Table 4(b). ^ the correspondence of each data set to a model can be readily observed. 

Although for each model there is an improvement with an increase in the number of coefficients, 
the difference is generally smaller for those groups of observations with higher average residuals. 

WEIGHTING OF IGRF 

It was decided that a weighted average of coefficients would provide the best compromise for an 
IGRF. Due to the restriction that models to be included should be based on truly spherical coeffi- 
cients, the RGO(LME) and USC&GS models were eliminated from the main field averaging. Since the 
surface-data residuals were so greatly influenced by crustal anomalies, it was decided to base these 
weights on the residuals to the satellite data. 

Several different weighting schemes were tried. Generally, the precise choice of weights used did 
not alter the overall results appreciably as long as those models best fitting the satellite data were given 
preference. The POG(X3/68) and AFCRL(1 1/67) models were eliminated from the considerations 
since each organization submitted another model. 



After the presentation of several semiquahtative arguments that the IGRF would be most useful 
circa 1965.0, the following tabulation of relative weiglits was agreed upon. Each weight was applied as 
an inverse square factor in combining the main field terms. 
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Model 


a (weight, in 7 ) 


GSFC( 12 / 66 ) 

40 

AFCRL(3/68) 

70 

RGO(3/68)-2 

80 

!ZMlRAN(3/68) 

100 


The GSFC model was given the 4 O 7 weight (even though it had a 6 I 7 residual to the (Xj 04 data) 
because the OGO 2 figure was 397 , the other satellite residuals were low, and it has the overall lowest 
residuals to the surface data. The AFCRL model and RGO contributions were roughly equivalent but 
the AFCRL was given a slightly smaller weight because of its lower residual to OGO 2, Cosmos 49, and 
the surface data. The IZMIRAN model was assigned a slightly higher weight because of its uncertainty 
in the polar regions, the model being derived from data at less than 50° latitude. This decision is sup- 
ported by the model’s relatively high residuals to data sets containing polar contributions (e.g., OGO 2, 
(XiO 4, observatory, land/sea, Scandinavian airborne, and Project MAGNET). 

There was less basis for rational comparisons in combining the secular change terms. Hence each 
model previously used was weighted equally, and the USC&GS and RGOl models were included since 
the secular change was independently derived for each. 

Although more lengthy considerations may have resulted in an improved procedure for deriving 
the first IGRF, this formulation provided a model composed of some contribution from each organiza- 
tion. At the same time, within the restrictions on the number of coefficients, it produced a model 
which agrees tolerably well with the test data set. This agreement is seen in the last column of Table 4. 
Surprisingly, the procedure appeared to produce a residual equal to or less than that of the contribut- 
ing models for some of the data sets at the n* =S truncation level. 

THE RESULTING MODEL IGRF(10/68) 

Since the IGRF is a composite of several models, it can be compared with each, as in Table 5. 
Here is listed for each of the contributing coefficient sets the deviation from the resulting IGRF. 
Although the disagreements between the various terms are sometimes relatively large for those with 
amplitudes of the order of 1 7 to lOy, those of higher magnitude are surprisingly close to one another. 
Of the main field terms, the largest discrepancy seems to be am^^ng those having m = n . 

The final IGRF( 10/68) coefficients are given in Table 6 . Maps of the field and its secular change 
are given in Appendix 1 . Appendix 2 gives a possible minor modification based on a suggested change 
of scale to a standard mean earth radius of 6371 in place of 6371.2 km. 

OBSERVATIONS AND RECOMMENDATIONS 

It is appropriate at this point to make some observations on the domain of applicability of the 
international reference field and on its deficiencies and limitations. As can be seen' in this report 
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Max. value of n 














h 


n 

m 

GSFC ( 12 / 66 ) 

AFCRL ( 3 / 68 ) 

RGO -2 

IZMIRAN 

USCGS 

RGO -1 

IGRF 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

16 

12 

-52 

11 

32 

5752 

2 

0 

0 

0 

0 

0 

0 

0 

0 

2 

1 

4 

-11 

6 

-9 

-6 

4 

-2006 

2 

2 

-11 

10 

-15 

71 

-7 

12 

130 

3 

0 

0 

0 

0 

0 

0 

0 

0 

3 

1 

G 

-6 

-22 

12 

-3 

-2 

-403 

3 

2 

-2 

-3 

-1 

22 

-3 

7 

242 

3 

3 

10 

-10 

20 

-53 

16 

-8 

-176 

4 

0 

0 

0 

0 

0 

0 

0 

0 

4 

1 

1 

-6 

12 

-10 

4 

-10 

149 

4 

2 

6 

-6 

0 

-28 

5 

-8 

-280 

4 

3 

3 

-7 

3 

-9 

7 

-2 

8 

4 

4 

-16 

-4 

11 

90 

35 

-7 

-265 

5 

0 

0 

0 

0 

0 

0 

0 

0 

5 

1 

-1 

11 

2 

-22 

4 

1 

16 

5 

2 

1 

3 

8 

-23 

1 

-5 

125 

5 

3 

-5 

5 

1 

23 

-4 

-7 

-123 

5 

4 

1 

-10 

10 

1 

9 

-3 

-107 

5 

5 

4 

3 

-2 

-25 


-7 

77 

6 

0 

0 

0 

0 

0 

0 

0 

0 

6 

1 

2 

-2 

-3 

-4 

3 

9 

-14 

6 

2 

-1 

-1 

1 

* 

6 

1 

3 

106 

6 

3 

1 

0 

-10 

8 

0 

0 

68 

6 

4 

2 

6 

1 

-26 

-12 

-2 

-32 

6 

5 

-2 

3 

-6 

15 

13 

-4 

-10 

6 

6 

3 

-4 

3 

-17 

-12 

-4 

-13 

7 

0 

0 

0 

0 

0 

0 

0 

0 

7 

1 

-2 

5 

11 

-16 

-6 

0 

-57 

7 

2 

0 

-3 

2 

0 

1 

2 

-27 

7 

3 

2 

-2 


-6 

-1 

-8 

-8 

7 

4 

2 

-5 

-1 

3 

-2 

-4 

9 

7 

5 

-1 

-3 

1 

8 

1 

-3 

23 

7 

6 

-1 

3 

-3 

3 

0 

-5 

-19 

7 

7 

0 

1 

-4 

4 

-2 

-5 

-17 

8 

0 

0 

0 

0 

0 

0 

0 

0 

8 

1 

-1 

7 

-7 

1 

7 

6 

3 

8 

2 

0 

-1 

6 

-9 

0 

7 

-13 

8 

3 

-1 

7 

-4 

-3 

2 

13 

5 

8 

4 

-2 

3 

1 

6 

5 

-4 

-17 

8 

5 

1 

-2 

1 

-6 

-5 

1 

4 

8 

6 

0 

-3 

2 

4 

4 

-7 

22 

8 

7 

0 

-1 

7 

-7 

-7 

-4 

-3 

8 

8 

-3 

3 

1 

8 

6 

-9 

-16 


TABLE 5(b) 


The IGRF w as developed to fill the need for a standard field model in which the permanence of a 
standard over a period of years outweighs the advantages of a high accuracy. Thus, the ultimate use of 
this model and fuiiher requests for revisions must be left to the users. 

The way to test the applicability of 1GRF( 10/68) to a particular problem is to perform regular 
tests of newer or more accurate models and compare the results with those based on IGRF. As the 
core field deviates more and more from the IGRF estimate, the accuracy will continuously decrease. 
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n 

m 

GSFC ( 12 / 66 ) 

AFCRLO / 68 ) 

RGO -2 

IZMIRAN 

USCGS 

RGO -1 

IGRF 

1 

0 

- 1.5 

- 3.3 

1.7 

3.3 

- 0.3 

0.2 

15.3 

1 

1 

0.7 

- 1.9 

0.3 

- 1.0 

1.8 

- 0.4 

8.7 

2 

0 

0.3 

- 3.3 

0.4 

- 0.7 

3.3 

- 2.2 

- 24.4 

2 

1 

- 0.5 

- 0.5 

- 1.3 

- 0.2 

2.4 

- 1.6 

0.3 

2 

n 

- 4.2 

- 2.0 

3.6 

0.4 

2.1 

2.9 

- 1.6 

3 

0 

- 1.7 

0.1 

- 0.2 

0.6 

1.1 

- 0.2 

0.2 

3 

1 

0.0 

0.6 

0.8 

1.9 

- 3.3 

1.3 

- 10.8 

3 

2 

1.7 

0.5 

- 1.7 

- 0.7 

0.2 

- 2.6 

0.7 

3 

3 

- 3.1 

- 3.0 

2.8 

2.1 

1.3 

3.2 

- 3.8 

4 

0 

2.1 

1.4 

- 0.3 

0.9 

- 4.1 

1.3 

- 0.7 

4 

1 

0.5 

0.6 

- 0.2 - 

0.1 

- 1.0 

0.8 

0.2 

4 

2 

1.3 

2.6 

- 1.0 

0.7 

- 3.6 

0.8 

- 3.0 

4 

3 

0.8 

- 0.4 

1.1 

0.1 

- 1.7 

0.3 

- 0.1 

4 

4 

- 1.4 

1.2 

- 0.9 

0.8 

0.3 

- 0.9 

- 2.1 

5 

0 

- 0.0 

1.4 

- 0.9 

- 1.3 

0.7 

- 1.1 

1.9 

5 

1 

- 0.5 

- C .5 

- 1.1 

- 0.8 

2.9 

- 0.7 

1.1 

5 

2 

0.0 

- 0.2 

0.1 

- 1.8 

0.9 

- 1.3 

2.9 

5 

3 

- 0.6 

- 0.2 

- 1.6 

- 1.0 

3.2 

- 0.9 

0.6 

5 

4 

- 0.6 

1.5 

- 1.0 

- 0.7 

0.8 

- 1.1 

0.0 

5 

5 

0.8 

2.7 

- 0.3 

- 1.2 

- 2.1 

0.4 

1.3 

6 

0 

- 0.4 

rl.l 

0.1 

- 0.4 

1.8 

0.1 

- 0.1 

6 

1 

1.2 

- 0.7 

0.3 

0.6 

- 1.3 

- 0.2 

- 0.3 

6 

2 

- 0.2 

- 0.6 

- 0.1 

- 0.4 

1.2 

0.7 

1.1 

6 

3 

0.7 

- 0.7 

0.1 

- 0.6 

0.4 

- 0.5 

1.9 

6 

4 

1.2 

0.9 

0.4 

0.8 

- 3.4 

1.0 

- 0.4 

6 

5 

0.6 

- 2.3 

0.4 

0.4 

0.9 

0.6 

- 0.4 

6 

G 

0.6 

0.6 

- 0.8 

0.2 

- 0.6 

- 2.0 

- 0.2 

7 

0 

- 0.1 

0.7 

0.5 

0.5 

- 1.6 

0.5 

- 0.5 

7 

1 

- 0.0 

- O.C 

0.3 

0.3 

- 0.6 

0.3 

- 0.3 

7 

2 

- 0.9 

- 0.3 

0.7 

0.7 

- 0.2 

0.7 

- 0.7 

7 

3 

- 0.4 

1.3 

0.5 

0.5 

- 1.9 

0.5 

- 0.5 

7 

4 

- 0.3 

- 0.8 

- 0.3 

- 0.3 

1.5 

- 0.3 

0.3 

7 

5 

- 0.7 

- 0.7 

0.0 

0.0 

1.4 

0.0 

- 0.0 

7 

6 

0.0 

0.9 

0.2 

0.2 

- 1.3 

0.2 

- 0.2 

7 

7 

- 0.1 

- 1.6 

0.6 

0.6 

0.5 

0.6 

- 0,6 

8 

0 

0.3 

- 0.3 

- 0.1 

- 0.1 

0.3 

- 0.1 

0.1 

8 

1 

0.1 

- 0.1 

- 0.4 

- 0.4 

0.8 

- 0.4 

0.4 

8 

2 

1.3 

- 0.5 

- 0.6 

- 0.6 

0.3 

- 0.6 

0.6 

8 

3 

- 0.2 

- 0.8 

- 0.0 

- 0.0 

1.1 

- 0.0 

0.0 

8 

M 

0.5 

- 0.8 

0.0 

0.0 

0.2 

0.0 

- 0.0 

8 

5 

- 0.0 

- 0.1 

r.i 

0.1 

- 0.0 

0.1 

- 0.1 

8 

6 

0.2 

- 0.7 

- 0.3 

- 0.3 

1.1 

- 0.3 

0.3 

8 

7 

0.1 

- 1.4 

0.3 

0.3 

0.6 

0.3 

- 0.3 

8 

8 

0.0 

- 0.4 

0.5 

0.5 

- 0.6 

0.5 

- 0.5 


TABLE 5(c) 


We have already made this test for the application to analysis of the time variations of the 
Cosmos 49, (XjO 2, and (XK) 4 data. For such studies the IGRF is not useful, the GSFC(1 2/66) 
model is insufficient, and fits based on the data themselves are being used. For higher accuracy stud- 
ies, we suggest using the GSFC( 12/66) model over the range 1900-1965 and the POGO( 10/68) model 
for 1965-1968. Beyond 1968, PCXXX 10/68) can be used until it is updated by more recent data and 
planned improvements in the analysis. 
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M 

n 

m 

CSFC ( 12 / 66 ) 

AFCRL ( 3 / 68 ) 

RGO -2 

IZMIRAN 

USCGS 

RGO-l 

IGRF 

1 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 

1 

- 1.7 

1.4 

2.3 

- 0.2 

- 1.7 

2.9 

- 2.3 

2 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

1 

- 2.3 

- 4.1 

- 0.2 

2.0 

4.7 

0.4 

- 11.8 

2 

2 

0.0 

1.3 

- 1.3 

0.9 

- 1.1 

- 1.5 

- 16.7 

3 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

1 

1.5 

0.8 

0.8 

0.1 

- 3.1 

- 1.0 

4.2 

3 

2 

1.8 

1.4 

0.3 

0.0 

- 3.7 

0.9 

0.7 

3 

3 

1.1 

- 3.3 

- 1.3 

2.1 

1.5 

- 0.8 

- 7.7 

4 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

4 

1 

- 2.1 

- 1.4 

1.1 

1.5 

1.0 

3.1 

- 0.1 

4 

2 

- 1.4 

- 3.0 

- 0.6 

- 0.7 

5.7 

- 2.3 

1.6 

4 

3 

- 1.2 

- 0.2 

0.1 

- 0.4 

1.7 

- 0.2 

2.9 

4 

4 

- 2.6 

- 0.3 

1.2 

2.6 

- 0.8 

1.5 

- 4.2 

5 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

1 

- 0.3 

- 0.1 

- 0.3 

- 1.9 

2.6 

- 0.4 

2.3 

5 

2 

- 0.1 

0.5 

0.3 

- 0.3 

- 0.3 

0.6 

1.7 

5 

3 

- 0.3 

0.6 

0.4 

1.2 

- 1.9 

0.6 

- 2.4 

5 

4 

- 0.3 

- 2.6 

0.2 

- 0.6 

3.4 

0.6 

0.8 

5 

5 

0.1 

0.2 

0.3 

0.5 

- 1.1 

0.8 

- 0.3 

6 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

6 

1 

1.1 

0.5 

- 1.1 

0.5 

- 1.0 

- 1.3 

- 0.9 

6 

2 

0.4 

0.5 

1.4 

0.6 

- 2.9 

0.6 

- 0.4 

6 

3 

0.6 

- 0.1 

- 0.0 

- 1.1 

0.8 

- 1.1 

2.0 

6 

4 

- 0.3 

1.3 

0.1 

l.l 

- 2.1 

- 1.0 

- 1.1 

6 

5 

0.4 

0.5 

- 0.1 

- 0.1 

- 0.7 

- 0.2 

0.1 

6 

6 

- 0.1 

- 1.3 

0.1 

- 0.9 

2.3 

- 1.1 

0.9 

7 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 

1 

0.1 

0.5 

1.1 

1.1 

- 2.7 

1.1 

- 1.1 

7 

2 

- 0.2 

- 1.1 

- 0.3 

- 0.3 

1.9 

- 0.3 

0.3 

7 

3 

0.1 

0.5 

- 0.4 

- 0.4 

0.2 

- 0.4 

0.4 

7 

4 

0.6 

- O.C 

- 0.2 

- 0.2 

0.4 

- 0.2 

0.2 

7 

5 

- 0.8 

- 0.8 

- 0.4 

- 0.4 

2.4 

- 0.4 

0.4 

7 

C 

0.4 

1.5 

- 0.2 

- 0.2 

- 1.4 

- 0.2 

0.2 

7 

7 

0.7 

1.1 

- 0.3 

- 0.3 

- 1.2 

- 0.3 

0.3 

8 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

8 

1 

- 0.7 

- 0.4 

- 0.1 

- 0.1 

1.3 

- 0.1 

0.1 

8 

2 

- 0.1 

- 0.0 

0.2 

0.2 

- 0.2 

0.2 

- 0.2 

8 

3 

0.4 

- 0.7 

0.3 

0.3 

- 0.4 

0.3 

- 0.3 

8 

4 

- 0.7 

- 0.2 

0.2 

0.2 

0.5 

0.2 

- 0.2 

8 

5 

0.4 

- 0.1 

0.3 

0.3 

- 1.0 

0.3 

- 0.3 

8 

6 

0.5 

0.2 

0.4 

0.4 

- 1.5 

0.4 

- 0.4 

8 

7 

- 0.1 

- 1.0 

0.3 

0.3 

0.5 

0.3 

- 0.3 

8 

8 

- 0.2 

- 0.4 

0.3 

0.3 

0.1 

0.3 

- 0.3 


TABLE S ( d ) 


The magnetic field derived from the IGRF or other magnetic field coefficients can be calculated 
from a wide variety of computer programs currently available. One such set of programs, based on a 
code originally developed by Jensen and Whitaker ( 1 960), may be obtained from 
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World Data Center A for Rockets and Satellites 
Goddard Space Flight Center (601 ) 
Greenbelt, Maryland 2077 1 


Tabic 6-Final IGRF( 10/68) coefficients. 


PPOrn - 1965.0 I , O . P . P .( l 0 / 68 ) 


n 

m 

f! 

h 

g 

h 

1 

0 

-30339 

0 

15.3 

0.0 

1 

1 

-2123 

5758 

8.7 

- 2.3 

2 

0 

-1654 

0 

- 24.4 

0.0 

2 

1 

2994 

-2006 

0.3 

- 11.8 

2 

2 

1567 

130 

- 1.6 

- 16.7 

3 

C 

1297 

0 

0.2 

0.0 

3 

1 

-2036 

-403 

- 10.8 

4.2 

3 

2 

1289 

242 

0.7 

0.7 

3 

3 

843 

-176 

- 3.8 

- 7.7 

4 

0 

958 

0 

- 0.7 

0.0 

4 

1 

805 

149 

0.2 

- 0.1 

4 

2 

492 

-280 

- 3.0 

1.6 

4 

3 

-392 

8 

- 0.1 

2.9 

4 

4 

256 

-265 

- 2.1 

- 4.2 

5 

0 

-223 

0 

1.9 

0.0 

5 

1 

357 

16 

1.1 

2.3 

5 

2 

246 

125 

2.9 

1.7 

5 

3 

-26 

-123 

0.6 

- 2.4 

5 

4 

-161 

-107 

0.0 

0.8 

5 

5 

-51 

77 

1.3 

- 0.3 

6 

0 

47 

0 

- 0.1 

0.0 

6 

1 

60 

-14 

- 0.3 

- 0.9 

6 

2 

4 

106 

1.1 

- 0.4 

6 

3 

-229 

68 

1.9 

2.0 

6 

4 

3 

-32 

- 0.4 

- 1.1 

6 

5 

-4 

-10 

- 0.4 

0.1 

6 

6 

-112 

-13 

- 0.2 

0.9 

7 

0 

71 

0 

- 0.5 

o.u 

7 

1 

-54 

-57 

- 0.3 

- 1.1 

7 

2 

0 

-27 

- 0.7 

?.3 

7 

3 

12 

-8 

- 0.5 

0.4 

7 

4 

-25 

9 

0.3 

0.2 

7 

5 

-9 

23 

- 0.0 

0.4 

7 

6 

13 

-19 

- 0.2 

0.2 

7 

7 

-2 

-17 

- 0.6 

0.3 

8 

0 

10 

0 

0.1 

0.0 

8 

1 

9 

3 

0.4 

0.1 

8 

2 

-3 

-13 

0.6 

- 0.2 

8 

3 

-12 

5 

0.0 

- 0.3 

8 

4 

-4 

-17 

- 0.0 

- 0.2 

8 

5 

7 

4 

- 0.1 

- 0.3 

8 

6 

-5 

22 

0.3 

- 0.4 

8 

7 

12 

-3 

- 9.3 

- 0.3 

8 

8 

6 

-16 

- 0.5 

- 0.3 


These programs internally convert the Schmidt-normalized coefficients to a more efficient Gauss- 
normalized form, update them to the epoch requested, and compute the geocentric components from 
the scalar gradient of the potential function, given the geocentric position. Conversions are also pro- 
vided from geodetic position to geocentric, as well as routines for rotating the output geocentric com- 
ponents into geodetic directions. Ignoring the differences between geodetic and geocentric coordinates 
will create errors up to about 200y. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Grccnbelt, MaryUii«{. Ja'*uary 13. 1971 
84 M 2-04- 12-5 1 
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Appendix 1 


Main Field Component and Isoporic Charts 
Computed From IGRF( 10/68) for 1965.0 at the Earth's Surfare 


The following figures represent the surface contours of the various geodetic components of the 
geomagnetic field and its secular change as computed by the IGRF. These diagrams are very similar to 
those given by Cain and Hendricks (1968) for the GSFC( 12/66) field and are drawn automatically 
using a computer program originally used for weather maps (Cain and Neilon, 1963). 

The plots are thus drawn to include the algebraic “lows” and “highs” of the component being 
displayed. These extrema occur at the center of the “+” or symbols. The dip poles arc noted for 
the H chart as “®”. 
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SECULAR CHANGE OF HORIZONTAL FORCE (gamma/year) 
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Appendix 2 


Coefficient Normalization 


All of the previous field derivations have arbitrarily set the earth’s mean radius at 6371 .2 for the 
value offl in the factors of the potential expansion. This value stemmed from the old standard 

earth constants with equatorial radius 6378.388 and flattening 1/297. However, the new constants 
have become 6378.165 and 1/298.25 respectively. Integrating 


we obtain 


where 




r cos Q dd , 





a - equatorial radius, 
and 

6 = d( 1 -/) is the polar radius with / the flattening factor. 

The values with the old and new constants are as follows. 

fab? 

297 6378.39 6356.91 6371.21 

298.25 6378.16 6356.77 6371.02 

It is recommended that for the sake of simplicity and not to be restricted to constants of only 
historical significance, we adopt the value of 637 1 for a. This is a very slight change and has only the 
effect of altering the term from -30,339 to -30,342 and the h\ term from 5758 to 5759. The con- 
stants a used to make the correction g = g' + otg\ where g' is the old value of g or /z, are 

n 123456789 10 11 12 

aX 1()5 9 13 16 19 22 25 28 31 35 38 41 44 

PRECEDING PAGE BLANK NOT FILMED 
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